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Effects of randomness on non-integer power law tails in multiplicatively 
interacting stochastic processes are investigated theoretically. Generally, 
randomness causes decrease of the exponent of tails and the growth 
rate of processes. Explicit calculations are performed for two examples: 
uniformly distributed and two peaked systems. Significant influence is 
demonstrated when a bare growth rate is low and coupling is weak. 
It should be emphasized that even the sign of the growth rate can be 
changed from positive to negative growth. 
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Power laws are widely observed not only in natural systems but also in social phe- 
nomena pin]. A number of efforts have been made to explain these behavior IBllHj but 
the underlying physics has not been clarified yet. Recently, inelastic Maxwell models 
where randomly chosen particles collide inelastically with each other have been studied 
extensively in the context of granular materials pillj. Then it has been shown that 
probability distribution functions of particle velocities have a non-integer power law 
tail asymptotically. The exponent of tails is determined from a transcendental equation 
and, generally, is a continuously varying function of a parameter, i.e., a coefficient of 
restitution. The system exhibits a kind of self-organized criticality. Non-integer power 
laws emerge not at a special (critical) value of a parameter but in a wide range of 
a parameter. ben-Avraham et al. have extended these studies to multiplicatively 
interacting stochastic processes. Here the system is composed of A^(^ 1) particles with 
positive variables Xi > 0{i = 1,2, ■ ■ ■ ,N) and evolves with multiplicative interactions 
between randomly chosen particles i and j 
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where a and b are positive parameters. When the eigenvalue Aa = a + & of the matrix 
A is more (less) than unity, the system grows (declines) totally. The case a + b = 1 
corresponds to inelastic Maxwell models where the total amount of variables (momenta) 
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is conserved. Power law tails appear ata + 6< 1, a > 1, or6> 1 and the exponent 
varies continuously with a and b. 

In actual processes, interactions sometimes accompany randomness. Especially, 
mesoscopic and/or macroscopic systems such as granular meterials are composed of 
not identical but polydisperse particles. In social phenomena, moreover, individuality 
of agents inevitably results in randomness in interactions. Note that the process is 
a reasonable model describing wealth distribution in economical systems ^21- In this 
Letter, therefore, we examine effects of randomness on power law tails in multiplica- 
tively interacting systems. 

We consider processes where a and b are random parameters with probability dis- 
tribution p{a,b). A master equation for probability distribution functions f{x,t) is 
given by 



^ + f = J J dadbp{a, b) 

dx'dy'-f{x',t)f{y',t){S{x — ax' — by')+S{x — bx' — ay')} . (2) 



The Fourier transform of eq.® is 
dg{k,t) 



dt 



+ g{k,t) = J J dadbp{a,b)g{ak,t)g{bk,t), (3) 



where g{k,t) = J dxf{x,t)e~^^^. The important property of eq.Q is that its moment 
equations become a closed set piTT]. Expanding as g{k,t) = 1 + J2'^=ifJ'n(t){—ik)"'/n\ 
and substituting into eq.Q, we have 



- KfJ'n = J J dadbp{a, b) ^ ^\(^^'_ ^-^ , /^^/^n-^Q V 



\n= f /dad6p(a,6)(a" + 6" - 1) =^ + F- 1. (5) 



Here the bar denotes average with respect to p{a,b). Equation (jH) can be solved 
sequentially, which enable us to treat eq.Q analytically. Here we pursue a similarity 
(scaling) solution of the type 



g{k,t) = <l>{v), {V = ke^') (6) 
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where 7 expresses a growth rate of the system. Substitution of eq.® into eq.© leads 
to 



d$(?7) 
IV- 



$(77) = J J dadbp{a,b)^{ar])^{br]). (7) 



di] 

The inverse Fourier transform of gives scaled distribution functions 



vl>(0 = /(x,t)e^* (e = xe- 



■■yt\ 



We pay our attention to tails in the small 77 and large ^ limit. This limit is classified 
into two cases and treated separately. In case I, we make the following ansatz, 



<^=l + Y.-^i-iriY + C^i-iri)% {l<n<a<n + l) (9) 

where a represents a non-integer power. Substituting eq.Q into eq.© and equating 
terms of each powers r]^, we get 



0{v) ■■■ 7 = Ai, (10) 



i-i 



^-^m!(£-m)! 

(£ = 2,3,---,n) 

Oiv'^) ■■■ ^a = X,. (12) 



The growth rate 7 is given by eq.()10p and the exponent a is determined from the 
combined equation of eqs.lfTIH) and (fT^ 



{a + h-l)s = a' + ¥-l. (13) 



The transcendental equation (fTSj) always has a trivial solution s = 1. When a nontrivial 
solution s = s* exists in a region s > 1, it gives the exponent a = s*. In case II where 
s* < 1, we make the ansatz, 



$ = l + a(-^^)^ (^<i) (14) 



3 



instead of eq.(jni). In this case, we have only the 0{'r]'^) term, eq. lfT^ . and the transcen- 
dental equation 



7S = a'* + 6* - 1. (15) 



Here we need one more condition to determine two variables 7 and s from one equation 
(I15|) . The answer is the selection of the minimum growth rate originally reported by 
ben-Avraham et al. [T^]. When the line 7s and the curve liF +¥ — 1 come in contact 
at s = s**, the growth rate 7 becomes minimum and s** gives the exponent a less than 
unity. 

Now we study effects of randomness. First, we compare the exponent a and the 
growth rate 7 with those cxo and 70 in the absence of randomness, which are determined 
from 



[s > 1) (16) 



7o = ao + &o - 1, 
(ao + bo- l)s = ag + 6g - 1, 

70S = ag + 6^ - 1, (s < 1) (17) 



where ao = a and bo = b. When s > 1, a* is a downwards convex function of s and 
a'^ > a**, while a* is upward convex and a* < a^ at s < 1. It follows that solutions of 
eqs. lfTSj) and (fT^ are always smaller than those of eqs. (fTB|) and (fTTj) . 



a < (To. (18) 



It becomes evident that randomess causes decrease of the exponent of tails. At the 
same time, the selection of the minimum growth rate leads to possible decrease of the 
growth rate 



7 < 70- (a < 1) (19) 



Next, we carry out explicit calculations for two examples. First example is an 
uniformly distributed system described by 



(ai < a < a2, bi < b < 62) 



p{a,b) = { («2-ai)(&2-&i)' • ■ (20) 

0. (otherwise) 
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Figure 1: Calculated values of (a) exponents a of tails and (b) growth rates 7 of 

processes in uniformly distributed systems at ao = 1.01, 60 = 0.001 ( ), ao = 1.01, 

60 = 0.01 ( ) and ao = 1.1, 60 = 0.001 ( ). 



The transcendental equations ()T3j) and ()T5j) read 



(ao + &o - l)-s 



7s 



1 



s + 1 
1 



, a2 - ai 
a2 — di 



+ 



62 - bi 



"2 



bl 



&2 - bl 



[s > 1) (21) 
[s < 1) (22) 



where ao = (ai+a2)/2 and bo = (6i+62)/2. Computed values of the exponent a and the 
growth rate 7 are plotted in Fig|T] as a function of A = (a2 — ai)/(2ao) = (62 — &i)/(26o)- 
When ao is close to unity and a bare growth rate of systems is quite low, a diminishes 
drastically with increasing randomness A, whereas decrease of 7 is remarkable at small 
bo and weak coupling. Notice that 7 changes even its sign from positive to negative 
growth. 

Second example is a two peaked process expressed as 



p(a, b) = p6{a — ai)6{b — 62) + (1 — p)S{a — a2)S{b — 62)- 



(23) 



Equation (|23|) gives rise to 



(ao + bo- l)s = p{al + bl) + (1 - p)(a^ + 6^) - 1, (s > 1) 



7S = p{al + bl) + (1 - p)(a^ + 6^) - 1, {s < 1) 



(24) 
(25) 
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Figure 2: Calculated values of (a) exponents a of tails and (b) growth rates 7 of 

processes in two peaked systems at j9 = 0.9 ( ), p = 0.5 (— — — ) and p = 0.1 

( ) with flo = 1-01 and 60 = 0.01. 



where Oq = pai + (1 —p)a2 and 60 = V^i + (1 — p)^2- Figure El shows numerical results, 
where A = (1 — p)(a2 — ai)/ao = {1 — p){h2 — hi) jho. Similarly in the case of uniform 
distribution, randomness A exerts serious influence at low bare growth rates and weak 
coupling. In addition, we find that when {1 — p) ^ 1, a few amount of impurities 
(aliens) with a large parameter 02 extremely reduce values of a and 7. 

From viewpoints of econophysics, conclusions obtained in this work might be im- 
portant for several reasons. (1) As mentioned before, individuality and associated 
randomness necessarily exist. (2) Low bare growth rates and weak coupling are realis- 
tic conditions. (3) The decrease of the exponent a of tails gives enlargement of wealth 
differentials and has a significant economical meaning as well as the decrease of the 
growth rate 7. It would be interesting to analyze wealth distribution by taking into 
account individuality and randomness. 



References 

[1] H. E. Stanley: Introduction to Phase Transitions and Critical Phenomena (Ox- 
ford, Oxford, 1971). 

[2] B. B. Mandelbrot: The Fractal Geometry of Nature (Freeman, San Francisco, 
1982). 

[3] R. N. Mantegna and H. E. Stanley: An Introduction to Econophysics (Cambridge, 
Cambridge, 2000). 

[4] H. Takayasu (Ed.): Empirical Science of Financial Fluctuations (Springer, Berlin, 
2002). 

[5] R. Albert and A.-L. Barabasi: Rev. Mod. Phys. 74 (2002) 47-97. 



6 



[6] P. Bak, C. Tang and K. Wiesenfeld: Phys. Rev. Lett. 59 (1987) 381-384. 

[7] H. J. Jensen: Self- Organized Criticality (Cambridge, Cambridge, 1998). 

[8] D. Sornette: Critical Phenomena in Natural Sciences (Springer, Berlin, 2000). 

[9] A. V. Bobylev, J. A. Carrillo and I. M. Gamba: J. Stat. Phys. 98 (2000) 743-773 
[10] M. H. Ernst and R. Brito: J. Stat. Phys. 109 (2002) 407-432. 
[11] E. Ben-Naim and P. L. Krapivsky: Phys. Rev. E66 (2002) 011309. 
[12] D. ben-Avraham, E. Ben-Naim, K. Lindenberg and A. Rosas: |cond-mat / 0308175 
[13] F. Slanina: |coiid-mat/0311235 



7 



